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Numerical calculation of the temperature rise in 
ignition by hot wires 


By Henrik L. SELBERG! and C. H. JoHANSSON? 


With 6 figures in the text 


The thermal conditions for a rapid rise in temperature in an explosive medium 
can in special cases be determined with satisfactory accuracy by analytical ap- 
proximations to the solutions of the differential equations. Frank-Kamenetsky 


thus successfully estimated the critical temperature for explosion. He calculated 


the maximum wall temperature at which a stationary distribution of tempera- 
ture is possible assuming decomposition according to the Arrhenius law and pro- 
portionality between the rates of decomposition and chemical heat production, 
The instantaneous rise in temperature which is required within a small volume, 
to obtain an explosion cannot, however, be estimated in this manner. This case 
is much more complicated and the only method available seems to be numerical 
integration with the aid of computors. In this way Cook (2) has calculated the 


time of ignition when a surface of an explosive stick is brought into contact 


with a hot bath kept at a constant temperature. 

This article gives the results of a numerical calculation of the temperature 
development in the explosive medium around a wire, the temperature of which 
is suddenly raised by applying a certain amount of energy. The explosive medium 
is assumed to decompose according to the Arrhenius law, when the explosive 
heat is released. 

The calculation has been carried out with the “Besk” (Stockholm) computer. 
The writers wish to express their thanks to Mr. G. Ehrling, Ph. L., who under- 


took the coding. 
Symbols 


Fi the time, cec. 

r the distance to the axis of the wire, cm. 

a radius of the wire, cm. 

T the temperature, K. 

7, the temperature in the wire and charge when ¢<0, K. 
7, the temperature of the wire when ¢=0, K. 

A the frequency constant, gar; 

£ activating energy, cal/mol. 
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R_ the gas constant, cal/mol K. 
Q heat of explosion, cal/g. 

n _undecomposed fraction of the explosive medium. 

0: density of the wire, g/cm. 

ct specific heat of the wire, cal/g C. 

0 density of the explosive, g/cm’. 

c specific heat of the explosive, cal/g C. 

A thermal conductivity of the explosive, cal/em s C. 

U_ energy per unit length of wire necessary for explosion, cal/em. 


The temperature of the wire is assumed to be the same across the entire ? 
cross section. The error resulting from this simplification is in all probability’ 
slight and decreases with a diminishing wire radius. The terminal effects are also) 
ignored, i.e. the wire is assumed to be infinitely long. 

For the temperature around the wire we have the equations 
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with the boundary conditions 


T=T) for r=oo 

T=T) and n=1 for t=0 and r=a 
T=T, for t=0 and r<a 

Peo = Loo for t>0 
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The calculation has been carried out for platinum wires 0.015, 0.030, 0.060 
and 0.120 mm in diameter surrounded by a composition with the same decom- 
posing properties as mercury fulminate. The following numerical values have 
been used: 


To = 293°K, 0;= 21.4 g/cm’, c,=0.033 cal/g C, 0 =2.5 g/em3, c=0.1 cal/g C, A= 
8.10 ° cal/ems C, A = 10" 5-1, H = 25400 cal/mol, Q@ =400 cal/g. The values of A 
and F are, according to Vaughan and Phillips (3). 

Whether the charge explodes or the heat supplied is led off from the wire 
before any considerable decomposition takes place must be decided from the tem- 
perature curves obtained. 

Fig. 1 shows the temperature versus time in the explosive for a = 15 pe (0.015 mm), 
r= 15.6 u and six equidistant initial temperatures between 7’, = 420 and 470°C. A 
rapid rise in temperature is obtained for 7',=470°C after 0.7-0.8 milliseconds. 
Fig. 2a shows the temperature versus time for 7, =470°C at 7 points in the ex- 
plosive medium between r=15.6 uw and r=19.2 4 at an interval of 0.6 pu. The ex- 
plosion occurs after 0.8-0.9 millisec. Fig. 2b reproduces the temperature according 
to fig. 2a as a function of the distance from the wire surface with the time as 
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Fig. 1. Temperature versus time 0.6 “4 from the surface of a wire 30 yu in diameter. The figures 
indicate the initial temperature of the wire. 
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Fig. 2a. Fig. 2b. 


Fig. 2a. Temperature versus time around a wire of 30 “ diameter and 470°C initial tempera- 
ture. The figures indicate rising radial distances with an interval of 0.6 uw. 
Fig. 2b. Temperature versus the distance from the surface of the wire according to fig. 2 a. 


parameter. Fig. 3 shows the undecomposed quantity as a function of the time un- 
der conditions according to fig. 2a. 

It is interesting to note that the temperature close to the wire surface after 
the first rapid rise can be practically constant for a long time before the tempera- 
ture rise leading to explosion begins. Fig. 4 gives an example of the develop- 
ment of the temperature showing a distinctive “induction period”. It is obvious 
that this is not a sign of any special delay in decomposition but is due to 
cooperation between the decomposition according to the Arrhenius function, the 
heat capacity and the heat transmission. 
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Fig. 3. Undecomposed quantity versus time on conditions according to fig. 2a. 
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Fig. 4. Temperature versus time 0.4 “ from the surface of a wire of 30 “ diameter and 470°C 
initial temperature. 


If 7. signifies the lowest wire temperature 7, at which an explosion occurs 
then 7’, is a function of the wire radius a. The following connection has been 
obtained. 


amm | Te° K | Oe = Te — 273 
0.06 683 410 
0.03 703 430 
0.015 743 470 
0.0075 823 550 


6.=T,.—273 is dotted in degrees Centigrade as a function of 1 /a in fig. 5. The 
diagram shows that the dots lie on the straight line 0. =390+0.12/a. According 


to this the explosion occurs when the energy supplied to the wire per em length 
exceeds 


U=1.11a+3437a% Ws/em (1) 


The continuous line in fig. 6 reproduces U as a function of a2 according to 
equation 1. 


The conditions for ignition by hot wires have been investigated by Jones (4). 
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Fig. 5. Fig. 6. 


Fig. 5. 0,=7',— 273 as a function of the inverted wire radius. 


Fig. 6. Minimum energy per em length as a function of a?. The unbroken curve represents eq. 1. 
The circlets are Jones’ measuring points for monoresorcinate. The dotted straight line represents 
eq. 2 and the broken curve U =3.81 a+ 3190 a? Ws/cm. 


The minimum energy per cm length of wire was measured at various strengths 
of current and extrapolated to instantaneous heating for a series of nichrome 
wires with diameters from 1.26 to 4.22 w. The explosive medium consisted of 
monoresorcinate, 80% by weight and chlorate of potassium, 20% by weight. 
Jones measuring-points are shown in fig. 6 as circlets. He assumed that they lay 
on a straight line with the equation 


U =0.0022 + 463002 Ws/em (2) 


and drew the conclusion that there exists a positive minimum energy charac- 
teristic of the explosive, below which ignition is impossible. The accuracy of the 
measurements is not sufficient, however, to decide if the functional connection 
for U is of the type in equation 1 or 2. It seems to us more probable that U 
tends to zero with the radius. 
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Tryckt den 15 april 1958 


Uppsala 1958. Almqvist & Wiksells Boktryckeri AB 
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